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Abstract 

In a large class of SUSY GUT models with see-saw mechanism of neutrino mass 
generation, lepton flavor violating (LFV) decays + + etc., are 

predicted with rates that are within the reach of present and planned experiments. 
A crucial element in these predictions is the matrix of neutrino Yukawa couplings 
Y(, which can be expressed in terms of the light and RH heavy neutrino masses, the 
neutrino mixing PMNS matrix U, and an orthogonal matrix R. Leptogenesis can take 
place only if R is complex. Considering the case of quasi-degenerate neutrinos and 
assuming that R is complex, we derive simple analytical expressions for the /i — > e + 7, 
r — ► ^+7 and r — ► e+7 decay rates. Taking into account the leptogenesis constraints on 
the relevant parameters we show that the predicted rates of the LFV decays \x — ► e + 7, 
and r — > e+7 are generically enhanced by a factor of ~ 10 3 to ~ 10 6 with respect to the 
rates calculated for real R, while the r — > \i + 7 decay rate is enhanced approximately 
by two orders of magnitude. 
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1 Introduction 



The solar neutrino experiments PjJ 121 El HI- the data on atmospheric neutrinos obtained 
by the Super-Kamiokande collaboration [5 , and the results from the KamLAND reactor 
antineutrino experiment 6J, provide very strong evidences for mixing and oscillations [311113 
of flavour neutrinos. The evidences for solar v e oscillations into active neutrinos z/ MjT , in 
particular, were spectacularly reinforced by the combined Super-Kamiokande and first SNO 
[3] data, by the more recent SNO data jlj, and by the just published first results of the 
KamLAND 6j experiment. 

The interpretation of the solar and atmospheric neutrino, and of the KamLAND data in 
terms of neutrino oscillations requires the existence of 3-neutrino mixing in the weak charged 
lepton current (see, e.g., [TUllTTj ): 



Here v\ L , I = e,/i,r, are the three left-handed flavor neutrino fields, Vj L is the left-handed 
field of the neutrino Uj having a mass rrij and U is the Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) neutrino mixing matrix [7 . It follows from the results of the 3 H /3-decay experiments 
[12J that rrij < 2.2 eV. The existence of the flavour neutrino mixing, eq. ((H), implies that 
the individual lepton charges, L e , and L T are not conserved (see, e.g., (EH)- Therefore, 
lepton flavour violating (LFV) processes like // — > e + 7, — > e~ + e + + e~, r — > ji + 7, 
fi~ + (A, Z) — > e~ + (A,Z), are allowed. However, if the neutrino (lepton) mixing in the 
weak charged lepton current is the only source of L e , and L T non- conservation, as in the 
minimally extended Standard Theory with massive neutrinos, the rates and cross-sections 
of the LFV processes are suppressed by the factor [T3j {mj/M w ) 4 < 5.6 x 1CT 43 , M w being 
the W ± mass, which renders them unobservable. 

The experimentally suggested smallness of the neutrino masses can naturally be ex- 
plained by the see-saw mechanism of neutrino mass generation ^3] . The see-saw mechanism 
requires the existence of heavy right-handed (RH) Majorana neutrinos. Right-handed neu- 
trinos (HI IE] are completely neutral under the Standard Theory gauge symmetry group. 
Consequently, they can acquire Majorana masses Mr that are not related to the electroweak 
symmetry breaking mechanism, and can, in principle, be much heavier than any of the known 
particles. The heavy RH Majorana neutrinos can generate through their CP-violating de- 
cays the observed baryon asymmetry of the Universe ^7j- I n grand unified theories (GUT) 
their masses are typically by a few to several orders of magnitude smaller than the scale of 
unification of the electroweak and strong interactions, Mx ~ 2 x 10 16 GeV. However, their 
presence in a theory can lead to a severe hierarchy problem associated with the existence 
of two very different mass (energy) scales: the electroweak symmetry breaking and the RH 
Majorana mass scale. In supersymmetric (SUSY) GUT theories the hierarchy between these 
two mass scales is stabilized. Hence, the SUSY GUT theories incorporating the see-saw 
mechanism of neutrino mass generation provide a consistent and appealing framework to 
account for neutrino masses and for the baryon asymmetry in the Universe. 

SUSY theories have additional sources of lepton charge non-conservation. In spite of the 
possible flavor-blindness of SUSY breaking, the supersymmetrization of the see-saw mech- 
anism, for instance, can induce new LFV effects ^H]- If SUSY is broken above the RH 



3 




(1) 



2 



Majorana mass scale, as, e.g, in gravity-mediated breaking scenarios, there are renormaliza- 
tion group effects that generate new lepton charge non-conserving couplings at low energy 
even if such couplings are absent at the GUT scale. In contrast to the non-supersymmetric 
case, these couplings give contributions to the amplitudes of the LFV decays and reactions 
which are not suppressed by the small values of neutrino masses and the LVF processes can 
proceed with rates and cross-sections which are within the sensitivity of presently operating 
and proposed experiments (see, e.g., [TH]). 

The solar and atmospheric neutrino data and the data from the reactor v e experiments 
KamLAND, CHOOZ and Palo Verde, were used successfully for determining the pattern of 
the 3 — v mixing and the values of the two independent neutrino mass-squared differences, 
Arrig and Am\ , which drive the solar and atmospheric neutrino oscillations. Under the 
rather plausible assumption of CPT-invariance, for instance, the recent KamLAND results 
practically establish jH] the large mixing angle (LMA) MSW solution as unique solution of 
the solar neutrino problem, with Am^ ~ 7 x 10~ 5 eV 2 and tan 2 9 Q ~ 0.40 favored by the 
data, #0 being the mixing angle which controls the solar v e oscillations. The analyses of the 
atmospheric neutrino data show that Am\ and the mixing parameter sin 2 29 a, responsible 
for the dominant atmospheric z/ M (z/ M ) oscillations into v T (P T ), have values Am\ ~ 3 x 
10~ 3 eV 2 ^> Am| , and sin 2 29a ~ (0.9 — 1.0). The existing data, however, does not allow 
to determine the sign of Am A . Furthermore, the neutrino oscillations are not sensitive to 
the absolute values of neutrino masses. Correspondingly, there are three different types of 
3-neutrino mass spectra which are compatible with the existing neutrino oscillation data [20] 
(see also, e.g., [2I]) : normal hierarchical (NH), mi <C m 2 -C m 3 , inverted hierarchical (IH), 
mi -C m 2 = m3, and quasi-degenerate (QD), mi ~ m 2 — m%, m\ 2 3 ^> Am\ . 

In the case of QD spectrum, neutrino masses can be measured directly in the 3 H (3- 
decay experiments which are sensitive to the v e mass, m Pe = m 12 ,3. The present bound 
obtained in these experiments reads [12], mi j23 = m Pe < 2.2 eV. Sensitivity to values of 
mi j2 .3 ~ 0.35 eV are planned to be reached in the KATRIN experiment [22]- If the massive 
neutrinos Vj are Majorana particles, as is predicted by the see-saw mechanism, neutrinoless 
double-beta decay experiments can also provide information on the type of the neutrino 
mass spectrum and on the absolute neutrino mass scale (see, e.g., [23 121] and the references 
quoted therein). They measure a combination of masses and mixing parameters known as 
the effective Majorana mass parameter, |<m>| (see, e.g., f3 E]). The most stringent 
constraints on |<m>| were obtained in the 76 Ge experiments: |<m>| < 0.35 eV [21] (90% 
C.L.), and |<m>| < (0.33 - 1.35) eV [25] (90% C.L.). Higher, or considerably higher, 
sensitivities to |<m>| are planned to be achieved in several (/3/3) 0l /-decay experiments of 
the next generation (for a review see, e.g., [2H])- H neutrinos have a QD mass spectrum, they 
can be relevant cosmologically through their contribution to the hot dark matter component 
of the Universe. The sum of neutrino masses (mi + m 2 + 777,3) can be determined with a 
precision of ~ (0.04 — 0.10) eV from cosmological and astrophysical data [2*7j . 

The Universe seems to be made only of matter; cosmologically significant amounts of 
antimatter have never been observed. This asymmetry between matter and antimatter can 
be understood as the result of the dynamical evolution of an initially symmetric Universe in 
which baryon number is not conserved, C- and CP- symmetries are violated and a devia- 
tion from thermal equilibrium exists |28j . If these conditions are fulfilled, baryogenesis, the 
process which generates an excess of baryons over antibaryons, can take place. At present, 
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one of the most favored scenarios for baryogenesis is the leptogenesis scenario jTTj in which 
the heavy RH neutrinos play a fundamental role. Their CP-violating and out-of-equilibrium 
decays produce a lepton asymmetry that is partially converted into a baryon asymmetry 
through anomalous electroweak processes. Leptogenesis has the attractive feature of provid- 
ing a link between neutrino masses and the baryon asymmetry. 

In a large class of SUSY GUT models with see-saw mechanism of neutrino mass gener- 
ation and flavour-universal soft SUSY breaking at the GUT scale (see, e.g., j2HI30j), the 
LFV processes and leptogenesis are related: they both depend (although in different ways) 
on the matrix of neutrino Yukawa couplings Y v . The latter is one of the basic ingredients 
of the see-saw mechanism. The matrix Y v can be expressed in terms of the light neutrino 
and heavy RH neutrino masses, the neutrino mixing PMNS matrix U, and an orthogonal 
matrix R. Leptogenesis can take place only if R is complex. Working in the framework of 
the indicated class of theories and taking R to be complex, we derive in the present article 
simple analytical expressions for the /x ^ e + 7, r — > /x + 7 and r — > e + 7 decay rates in the 
case of quasi-degenerate neutrino mass spectrum. We use the model of leptogenesis of [3T] , 
in which the heavy RH neutrinos are produced non-thermally in inflaton decays, to obtain 
constraints on the parameters which determine the leading contribution in the LFV decay 
rates. Taking into account the leptogenesis constraints, we show that the rates of the LFV 
decays /x — >• e + 7, r ^ /x + 7 and r — ► e + 7, obtained for complex R, are generically strongly 
enhanced with respect to those calculated for real R. We present quantitative results for the 
enhancement factors for the indicated three LFV decays. 

Detailed predictions for the rates of the LFV processes in the class of SUSY GUT models 
with see-saw mechanism considered in our work were obtained, e.g., in refs. EDI E2 
EH1 El]- The case of quasi-degenerate neutrinos we analyze was discussed, in particular, in 

E21 EH] • However, the results in these articles were obtained for real matrix R. In |J3] 
the case of hierarchical neutrino mass spectrum was considered. The articles quoted in ref. 
30J contain rather comprehensive study of the LFV processes, including the case of complex 
R and the leptogenesis constraints, but for light neutrino mass spectra with normal and 
with inverted hierarchy. 

2 The neutrino Yukawa coupling 

The superpotential of the lepton sector in the MSSM with RH neutrinos is given by: 
W lepton = It T Y e L ■ H d + N L C T Y U L ■ H u - ^NfM R N c L , (2) 

where the family indices were suppressed. Here Lj, j = e, /x, r = 1, 2, 3, represent the chiral 
super-multiplets of the SU{2) L doublet lepton fields, l c L , j = e, /x, r = 1,2,3, is the super- 
multiplet of the SU{2) L singlet lepton field l c . L = CW R , where C is the charge conjugation 
matrix and l jR is the right-handed charged lepton field, Nj L is the super-multiplet of the 
SU(2) L singlet neutrino field N? L = CNj R , where N jR is the RH neutrino field, and H u 

and Hd are the super-multiplets of the two Higgs doublet fields H u and Hd carrying weak 
hypercharges — | and ~, respectively. In eq. (0), Y v is the 3x3 matrix of neutrino Yukawa 
couplings, Y e is the 3x3 matrix of the Yukawa couplings of the charged leptons, and M R 
is the Majorana mass matrix of the RH neutrinos N jR . We can always choose a basis in 
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which both Y e and M R are diagonal. We will work in that basis and will denote by Dm the 
corresponding diagonal RH neutrino mass matrix, Dm = diag(Mi, M 2 , M3). 

The see-saw mechanism generates a Majorana mass matrix for the left-handed neutrinos 
of the form: 

m„ = (Y u v u ) T D^(Y u v u ), (3) 

where v u is the vacuum expectation value of H u . The neutrino mass matrix is diagonal- 
ized by a single unitary matrix U according to 

D m = U T m u U = diag(m 1 ,m 2 ,m 3 ) , (4) 

where U is the PMNS matrix in the weak charged lepton current, eq. (J"fJ). 

It is convenient to choose rnj > 0, to number the massive neutrinos in such a way that 
mi < vtl2 < iris, and to work with Majorana neutrino fields Uj which satisfy the Majorana 
condition: C{pj) T — fj, j — 1, 2, 3. In this case the PMNS matrix U can be written as 

U = V ■diag(l,e ia ,e ip ), (5) 



where a and {3 are two Majorana CP- violating phases [o5]. For V one can use the standard 
parametrization 

(C13C12 C13S12 s 13 e~ lS \ 

-C23S12 - S23Si 3 ci 2 e l<5 C23C12 - s 2 3S 13 s 12 e t5 S23C13 , (6) 
iS iS I 

S23 s 12 — C23Si3Ci2e — S23C12 — C23Si3Si2e C23C13 J 

with the usual notations, s»j = sin^-, etc. If, for instance, Affig = Am^ (neutrino mass 
spectrum with normal hierarchy) and Am\ = Am^, one can identify 6u = 9®, 823 = #a? 
while ^13 is limited by the data from the CHOOZ and Palo Verde experiments [361 13*7] . 
sin 2 13 < 0.05. 

The matrix D m can be expressed as 

D m = U T Y T u v u D- M l Y u v u U = U T Y T v v u D K } l2 D h l /2 Y v v u U. (7) 

Following ref . j2H] , we define the complex matrix R: 

R = D k ] l2 Y u v u UD m ^ . (8) 

Given D M 2 , D^ 2 and U, the most general neutrino Yukawa coupling matrix reads 

Y„ = -D M 2 KD^ . (9) 

It follows from eq. (J7J that R is an orthogonal matrix, RR T = 1. In order for the 
leptogenesis scenario of baryon asymmetry generation to work, R must be complex and we 
will keep R complex throughout this study. As we will see, apart from being a necessary 
condition for leptogenesis, this leads also to drastically different predictions for the rates of 
the LFV processes like \i — >e + 7, r— >e + 7, r^/i + 7. 

The see-saw model contains 18 physical parameters - 6 phases and 12 moduli. These 
include, in the basis we work, the 3 (real) masses of the heavy RH Majorana neutrinos, 
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and 9 moduli and 6 phases of Y u (3 of the 9 phases in Y u can be eliminated through a 
rephasing of the LH charged lepton fields). At low energies it is convenient to parametrize 
the model by the 3 angles and 3 phases of the PMNS mixing matrix U, eqs. (jSJ) and ((HJ), the 
3 light neutrino masses, m lj2 ,3, and the 6 parameters - 3 moduli and 3 phases, of the complex 
orthogonal matrix R. The 3 additional real parameters of the model are the 3 heavy RH 
Majorana neutrino masses, contained in Dm- 

The 3 lepton mixing angles in the PMNS matrix U and the two neutrino mass squared 
differences, Amg and Am\ , can be measured with a relatively high precision in neutrino 
oscillation experiments. These experiments could also provide information on the Dirac 
CP-violating phase 5, whereas information on the two Majorana CP-violating phases, a 
and (3, can be obtained, in principle, in processes in which the Majorana nature of neutrinos 
manifests itself, such as (/?/?) 0l ,-decay, K~ — > 7r + +/i~-|-//~ decay, etc. (see, e.g., [SHUSH])- The 
measurement of the neutrino mixing parameters would be complete with the determination 
of the type of the neutrino mass spectrum and of the absolute neutrino mass scale. 

The probabilities of the LFV processes and the baryon asymmetry in leptogenesis depend 
on the see-saw parameters respectively via the quantities 

Y+Y„ = \uD£B! D M RD^ , (10) 

v 

and 

Im [(Y.Yt)^.] 2 = llm [(D M *R D m & D 1 ^} 2 , i ± j. (11) 

u 

Thus, the matrix R enters into both the expressions for the rates of the LFV processes and 
for the baryon asymmetry. 

We will consider in what follows the case of quasi-degenerate neutrino mass spectrum, 
m\ = m 2 = m 3 , mf 2 3 >> Am\ , Am^ . We can then write 

mx = m v , m 2 = m u + — — Ami , m 3 = m u + -^—Am 2 A , (12) 

lm v lvn v 

where m v is the neutrino mass determining the absolute neutrino mass scale which is not 
known, m v < 2.2 eV |12J. It is natural to assume that also Dm has quasi-degenerate 
eigenvalues, M 12 ,3 — Mr, D m — M R 1; otherwise, an exceptional fine-tuning between Y u 
and D M would be needed in order to obtain a QD spectrum for the light neutrinos. 
In the case of QD neutrino mass spectrum one has 

Y v = —M^m^Rdiag^l + Aml /(4m£), 1+Ami /(4m^))C/ t = —M^m^RU^ . (13) 

Hereafter corrections 0(Am 2 & /(2ml)) anc ^ 0(Am\ /(2m 2 ,)) will be neglected. 
The matrix R can be parametrized as 

R = e iA 0, (14) 

where A and O are real matrices. The orthogonality of R implies that O is orthogonal and 
A is antisymmetric. A different parametrization of R has been used in previous works (see, 
e.g., [29 1), but this one is particularly useful if the neutrino mass spectrum is of the QD type. 
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Up to corrections of the order of Am A /(2ml) an d /(^ m l): i.e., in the approximation 

of exact degeneracy of the three Majorana neutrinos ^1,2,3, the matrix O can be absorbed in 
the PMNS matrix U - the latter is defined up to a real orthogonal matrix and U and U O lead 
to the same physics Thus, up to relatively small corrections, O can effectively be taken 
to be the unit matrix, O = 1. This simplification is due to an additional 0(3) symmetry 
present in the lepton sector when the neutrino mass spectrum is exactly degenerate [40]. 

Thus, up to corrections of the order of Am A /(2ml) an d Amg /(2ml), the matrix R in 
the expression for Y u , eq. f|13j) . is effectively given by e* A . The matrix e lA can be explicitly 
calculated in terms of the three non-zero elements of A. If we write 



A = -a c , (15) 




then 



cosh r — 1 . r, . sinh r 



e* A = 1 - A 2 + i A , (16) 



where r = \J 0? + b 2 + c 2 . 

Our final expression for Y v in the case of QD neutrino mass spectrum is 

Y„9i — M R /2 m x J 2 e iA tf , (17) 

with e tA given by (jlfij) . 

If we take the mixing angles 6*12 = 9 Q and #23 = Oa as known and neglect sin 6*13 in U, 
both YlY u and Y„Y^ depend on 5 real parameters: Mr, m u , a, b, c; YlY u depends in 
addition on the phases a and (3. 

Yukawa couplings are expected to have moduli less than one, \Y U \ < 1. Taking a = b = 
c = k we get from the diagonal elements of Y v the condition 



cosh3A; < 



261GeV 1 



y/M R m v 2 

so that k < {1.4, 0.9, 0.3} for m u = 0.2 eV and M R = {10 10 , 10 12 , 10 14 } GeV, respectively. 

3 The processes — > £j + 7 

The existence of two Yukawa couplings in the lepton sector generally causes lepton 
flavour violation in a way analogous 14, to its quark sector counterpart in the Standard 
Theory. In the minimally extended Standard Theory with massive neutrinos and in the 
non-supersymmetric versions of the see-saw model, the decay rates and cross sections of the 
LFV processes are extremely suppressed: one has, for example, for the branching ratio of 
the [A — > e + 7 decay, BR(fi — * e + 7) < 10 -47 [HI HI]. Such small branching ratios are 
unobservable. The present experimental limit is 02] 

BR(fj, -> e + 7) < 1.2 x HT 11 . (19) 
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This bound is expected to be improved at least by a few orders of magnitude in the future. In 
an experiment under preparation at PSI [131 > f° r instance, it is planned to reach a sensitivity 
to 

BR(n -> e + 7) ~ 1(T 14 . (20) 

As we have seen, the rates of the LVF processes in the minimally extended Standard 
Theory with massive neutrinos are so strongly suppressed that these processes are not ob- 
servable in practice. In a SUSY theory the situation is very different because there is a new 
source of lepton flavor violation: the soft SUSY breaking Lagrangian, S£ so j t - The breaking 
of SUSY will, generally, cause lepton flavor violation. Indeed, off-diagonal elements in the 
neutrino Yukawa coupling can give rise to off-diagonal elements in the slepton mass matrix 
at low energies through renormalization group effects. 

The slepton sector of the soft supersymmetry breaking Lagrangian has the form 

- Sf^ = (mDylrJZj + (mj)ije* Ri e Rj + {va^ij^Rj 

+ (A^H&nLj + A^HJj^Lj + h.c) . (21) 

Lepton flavor violation can be generated by off-diagonal elements of the soft SUSY breaking 
parameters. The most conservative starting point for J? so ft is the assumption of universality 
at the GUT scale M x : 

( m lh = ( m l)ij = ( m lh = s v m o > 

™H d = ™< 2 Hu = m 2 , (22) 
A u = Y u a m , A e = Y e a m . 

Thus, at the unification scale, flavor is exactly conserved by Jf S oft- Nevertheless, soft SUSY 
breaking terms suffer from renormalization via Yukawa and gauge interactions. In this way, 
LFV in the Yukawa couplings will induce LFV in the slepton mass matrices at low energy 
even if the slepton masses are flavour-universal at high energy. 

The RGE for the left-handed slepton mass matrix is given by (see, e.g., (2111301) 



d , 0s d 



1 



+ 16^ [KYJY* + Y J Y „mi 



where the first term is the standard MSSM term which has no LFV, while the second one is 
the source of LFV. In the leading-log approximation and with universal boundary conditions, 
the off-diagonal elements of the left-handed slepton mass matrix at low energy are given by 



m 



2 ^ - ^(6 + 2a2K (Y tY^>g^. (23) 



This equation shows the connection between LFV in neutrino Yukawa couplings and LFV 
in slepton mass terms. 

Let us turn now to the lepton-flavor violating processes of the type l{ — > Ej + 7. The 
amplitude for this process has the general form 

T = e a £,m^a af5 q^A L P L + A R P R )ii , (24) 
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Figure 1 : Feynman diagrams giving the dominant contribution to the 1^ — > £j + 7 decay amplitude in the 
mass-insertion approximation. \A denotes charginos and neutralinos, and Li are the slepton doublets in the 
basis in which the gauge interactions and the charged-lepton Yukawa couplings are flavour-diagonal. The 
photon can be emitted from the chargino or the slepton lines. 



where q is the momentum of the photon, Pr(l) = + )7s)/2 and Al (A r ) is the coefficient 
of the amplitude when the decaying lepton is left-handed (right-handed). The corresponding 
branching ratio is 

12tt 2 

BR(£i - £j + 7) = ^(\A L \ 2 + \A R \ 2 ). (25) 

The terms |Al,.r| contain the contributions of the neutralino and the chargino loops (see 
Fig. 1). Explicit expressions for Al and An can be found in the literature jSj. In the mass 
insertion approximation, the diagrams contributing to £j — ► £j + 7 have the generic form 
shown in Fig. ^ and the branching can be estimated using the expression 



BRft _> i j + J ) ~ ^\A R \ 2 ~ ^M^tan 2 /?, (26) 
Gp Gp rrig 



where mj represents a scalar lepton mass. In the leading- log approximation, using (|23|). one 
finds [29, 30 j 



-> ^ + 7) 



a 3 



m%G 2 F 



3 + a 2 2 M 



m log ■ 



2 



kYtY^/tan 2 /?. (27) 



Therefore, the off-diagonal elements of Yj Y„ are the crucial quantities needed to estimate 
the branching ratios. 

Using the expression for Y u in eq. ([17]) we find that in the case of QD neutrino mass 
spectrum and in the approximation of negligible splitting between the neutrino masses, 

{YtYJv ~ \M R m u (Ue l2A U% . (28) 

For small values of a, 6, and c, and negligible s 13 we obtain 

(Ue i2A U^) 21 -2t [-a(c\ 2 e ia + s 2 l2 e- ia )c 23 - e ip s 23 (bc 12 + cs 12 e- ta )] (29) 
(£/e i2A f/t) 31 ~ 2 2 Hc\ 2 e ta + s^e"**)^ - e i/3 c 23 (6c 12 + cs 12 e- ia )] (30) 
(Ue i2A U% 2 ~ 2* [-2ias 12 ci2C 23 S23 sin a + (&s l2 - cci 2 e ia ) (s 2 3 e i/3 + c 2 3 e" i/3 )] (31) 

These elements control the /i — > e + 7, r— > e + 7 and r — >■ ^ + 7 decay rates, respectively. 
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Results for the LFV decay rates of interest in the case of real R have been obtained, e.g., 
in |2"9~1 [32J We will analyze next the differences in the predictions for the fi — > e + 7, 
r— >e + 7, r— >/i + 7 decay rates which appear when the matrix R is complex. We will 
denote the matrix Y£Y„ obtained by taking R real as YjY„|]R. One has [2TH 132] 

/ j. x Mr 

(YtY„| H ) H , = [UnUt 2 (m 2 - m x ) + U l3 U^(m 3 - mi)] (32) 



K 
Mr 

II 



Ami Ami , , 

I = = e and I — r,l = e, yU.(33) 



2m u 3 2m u 



This expression does not depend on R and has the special property that it depends on the 
masses only through the differences m-i — m\ and m 3 — m\ but not through their absolute 
values. Thus, in contrast to eq. (|28j). there is no contribution proportional to m v . The main 
contributions in the cases of \i — > e + 7 and r — > e + 7 decays are of order s\ 3 Am\ / (2m,,), 
or of order Am 2 & / (2m u ) if S13 is rather small; in the case of r — > /i + 7 decay it is of order 
Am\ /(4m„). Therefore, as long as \a\, \b\, \c\ ^ 10~ 3 , the fx — > e + 7, r— > e + 7 and 
r — > /i + 7 decay rates calculated using eq. (J2BJ), will typically be enhanced with respect 
to the rates calculated using eq. (jSHJ), because in the case of QD neutrino mass spectrum 
one has m v 3> Am\ /(4m u ), s\ 3 Am\ /(2m 1/ ),Am| /(2m„). The precise magnitude of this 
enhancement depends on the values of the parameters a, b and c, contained in A. 

Expressions (129)1 - (J3*T|) and ()33|) differ significantly in one more aspect: in contrast to 
YJY^|k, eq. (jHSJ), the quantity Y^Y^ calculated for complex R depends on the Majorana 
CP-violating phases a and (3 in the PMNS matrix 2 U. Thus, the observation of the LFV 
processes \i — > e + 7, r — > /i + 7, etc. could allow one to get information about these phases. 
Let us recall that determining or even constraining the Majorana CP-violating phases in the 
neutrino matrix U is a formidable problem (see, e.g., [2*T| I38|). 

In the numerical estimates which follow we take s±2 = sin# Q = 0.6, S23 = sin#A = l/v2, 
5 = 0, Am\ = 3 x 10~ 3 eV 2 , Am| = 7 x 10~ 5 eV 2 , m v = 0.3 eV and consider two different 
values for s 13: 0; 0.2. We have 

IfYtYM l 2 ~lfYtvn | 2 ~. M >' y f 6xl0- 6 if Sl3 = 0.2, 
KY^kj^l ~ |(Y,Y,| R ) 21 | _ __ x j L4 x 1Q _ 8 ifsi3 = 0Q (34) 

We will compare these results with the results we get for a complex matrix R. We set 
a = 7r/2 and (5 = 0. If we choose |a|, |6|, |c| ~ (9(10 _1 ), we always get a result that is much 
larger than ()34j) . We have, for instance, 

irvtv^ 12^ / 0.34, for (a, 6, c) = (0.2,-0.4,0.5), 

K^^il - w 4 x I 0.81, for (a, 6, c) = (0.4,0.3,0.2) ' {60) 

where we used eq. (|25|i with s 13 = (the results for s 13 = 0.2 are only slightly different). We 
see that the coefficient in the right-hand side of the above equation is always O(0.1 — 1.0). 
This result means that the branching ratio of the \x — > e + 7 decay will be enhanced with 
respect to the prediction based on eq. (j54"jl 3 approximately by a factor of 10 5 to 10 8 



2 This is valid also in the cases of neutrino mass spectra with normal and inverted hierarchy |30| . 
3 That if R is complex, BR(ii — > e + 7) could be enhanced with respect to the branching ratio predicted 
for real R, was noticed in |29j . 
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depending on the value of S13. Even if we take \a\, \b\, \c\ ~ C(10~ 2 ), there is still an 
enhancement of about three to six orders of magnitude. The same results are valid for the 
r — > e + 7 decay rate. 

The r — > /i + 7 decay rate is also enhanced, but the magnitude of the enhancement is 
smaller than in the case of the \x — > e + 7 decay rate: by a factor of ~ 10 4 for \a\, \b\, |c| ~ 
0(10^), and of ~ 10 2 for \a\, |6|, |c| ~ O(10" 2 ). This is due to the fact that the leading 
term in (Y^Y^| R ) 32 is not suppressed by s 13 . 

Detailed predictions for the rate of the /1 — > e + 7 decay, obtained for real R and for 
QD neutrino mass spectrum, can be found in [SHI E21 01j • They can be used, together with 
eqs. and (jHEJ), to estimate BR(fi — > e + 7) for complex R we have considered. The 

substantial enhancement we have found certainly makes the importance of the searches for 
this decay even more significant. 

Within the see-saw model, the neutrino Yukawa couplings, Y„, plays a major role in the 
generation of neutrino masses and in determining the rates of the LFV processes such as 
/U^e + 7, r^yU + 7 and r — > e + 7 decays. It plays a fundamental role also in leptogenesis. 

4 The Leptogenesis Constraints 

The convenient dimensionless number which characterizes the magnitude of the baryon 
asymmetry of the Universe is the ratio of the baryonic charge density, n^ — rig, to the entropy 
density, s. The presently observed baryon asymmetry is 

Y B = nB ~ ns = (0.1 - 1) x 10- 10 . (36) 
s 

The aim of baryogenesis is to explain this number in terms of processes and fundamental 
parameters of particle physics. In leptogenesis, the out of equilibrium decays of heavy RH 
neutrinos produce a lepton asymmetry which is reprocessed by sphaleron processes into a 
baryon asymmetry. If the light neutrinos are quasi-degenerate, m v ^0.1 eV, the out-of- 
equilibrium condition cannot be satisfied and the amount of produced lepton asymmetry is 
strongly suppressed (see, e.g., jSJHHl), unless the RH neutrinos are produced non-thermally. 
We shall consider leptogenesis via decays of RH neutrinos Ni which are produced through 
inflaton decays |3T] . 

At tree level the decay width of a heavy neutrino iVj is, 

T m = T(Ni ^H u + l) + T(Ni ^H c u + V s ) = ^-(Y^YDuM, . (37) 

If CP is not conserved by the neutrino Yukawa couplings, the interference between the 
tree and the one-loop diagram contributions to the Ni decay amplitudes results in a lepton 
number production. The lepton number asymmetry per decay of a RH neutrino is 

tl ~ F(Ni H u + l) + T(Ni H° + l c ) 

- -^7^7V:E Im [{( Y ^)^ 2 ] [f(M]/M?)+g(M]/M?)} . 

V f V )%% j-^i 

(38) 
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Here H u , /, and Aj denote scalar or fermionic components of the corresponding supermulti- 
plets, / is the contribution from the one-loop vertex correction |47| 14*%] HE] 



f(x) = yfx 



log 



X 



X 



(39) 



and g is the contribution from the one-loop self energy diagrams, which can be reliably 
calculated in perturbation theory if the condition 

| Mi - Mj\ > |r € - Tj\ (40) 

holds. One finds 

9{x) = ^ ■ (41) 
x — 1 

For quasi-degenerate RH neutrinos, x ~ 1 and g ^> /. 

The ratio of the lepton number density to the entropy density s produced by the 
inflaton decay is given by [3T] 

— = IY, e * BR ^ -> N i N i)~ » (42) 

where <p denotes the inflaton field, BR((p — ► A^A^) = Sr^ is the — > A^Aj decay branching 
ratio, is the mass of the inflaton and Tr is the reheating temperature after the inflation. 
We have assumed that Mr > Tr in order to prevent lepton-number violating processes from 
washing out the lepton asymmetry after the A^'s have decayed. Part of the lepton asymmetry 
is immediately converted into baryon asymmetry via the sphaleron effect, 



^ = C 7 ^, (43) 



with C ~ -0.35 in the MSSM. 
Using eq. Q we obtain 



Y„Yt = ^DtfRDn&Dtf . (44) 



Hence, in general, leptogenesis is independent of the mixing angles and phases contained 
in the PMNS matrix U. If R is real, Im(Y„YJ) = and leptogenesis cannot work. This 
is a model independent statement and it is the main reason we have to assume that R is 
complex. 

For quasi-degenerate neutrinos, eq. ()44|) can be further simplified, 



Y Y t ^ M ^ m y 2A 



, , ■ (45) 

v 

u 



It is well-known that if the RH neutrinos are completely degenerate, the generated lepton 
asymmetry is zero [4*9*] . Thus, one has to break the exact degeneracy in the heavy RH 
neutrino masses. We write 

M 2 = M 1 (l-e 2 ), M 3 = M x (l - £ 3 ), |e 3 |»N- (46) 
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\e 2 \ > 


1 


k 2 


-b 2 \ 


Mr 


47T 


10 14 GeV ' 


N > 


1 


k 2 


-a 2 \ 


M R 


47T 


10 14 GeV ' 


N > 


1 


|6 2 


-a 2 \ 




47T 


10 14 GeV ' 



The maximal values of \e 2 \ and (£3! which are naturally consistent with a low-energy quasi- 
degenerate neutrino mass spectrum are \e 2 \ — Amg/2m^ and ks| — Am\ /2m 2 ,. 
Conditions (}4T)j) for small r translate into 

Mr, 

(47) 
(48) 
(49) 

where we have used m v = 0.3 eV and v u = 174 GeV. Since \e 2 \ — 10 -4 and \es\ ~ 10~ 2 , for 
the extreme value Mr ~ 10 14 GeV these conditions are satisfied as long as |c| ^ 10~ 2 and 
\a\ < 10" 1 . For Mr ~ 10 10 GeV, eqs. (gZ| - (03J) lead to the constraints |a|, |6|, |c| < 1. 
We shall compute next the lepton number asymmetries, eq. (|38jl . From (|TH|) we get 

Im [(e i2A ) 12 (e i2A ) 12 ] = 2^ sinh2r(cosh2r - 1) 

= Im [(e l2A ) 23 (e j2A ) 23 ] 

= Im[(e i2A ) 31 (e* 2A ) 31 ] (50) 

and Im [(e l2A ) iJ (e j2A ) iJ ] = — Im [(e* 2A ) Ji (e* 2A ) Ji ] . Assuming that a, b and c are small, we can 
expand the hyperbolic functions in ()50j) to obtain 

\M R m v abc 

ei = e 2 ~ 5 (51) 

2M R m u abc 

e.3 5 52 

71 vl e 3 

The baryon asymmetry thus generated is 



n B 
s 



where we have neglected the €3 contribution. Hence, the empirical baryon asymmetry is 
obtained with a reheating temperature of Tr ~ 10 8 GeV for \abc/e 2 \ ~ 10 _1 and a natural 
choice of the remaining parameters. Since \e 2 \ ~ 10~ 4 , we have 

|o6c| ~ 10" 5 . (54) 

Larger values of \abc\ are possible if 2M^/m0 <^ 1. Higher reheating temperatures would be 
compatible with smaller values of the product \abc\, but would also lead to the cosmological 
gravitino problem |50j . 

The baryon asymmetry is proportional to the product abc and therefore none of the 
three parameters can be zero. Moreover, they cannot be exceedingly small, otherwise it 
would be impossible to reproduce the baryon asymmetry. Equation (|34*j) implies that at 
least one of them - \a\, \b\, or |c|, must be of order 10~ 2 or larger. This number fixes the 
possible enhancement of BR(fi —>■ e + 7): about four orders of magnitude for s 13 = 0.2 and 
six orders of magnitude for s 13 = 0. BR(r — > e + 7) is enhanced by similar factors, while 
BR(t —rfi + y) is enhanced approximately by two orders of magnitude. 
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5 Conclusions 



We have considered the fj, — ► e + 7, r — > + 7 and r — > e + 7 decay branching 
ratios in a class of SUSY GUT models with see-saw mechanism of neutrino mass generation. 
We have assumed that the orthogonal R matrix which was introduced in and which is 
related to the neutrino Yukawa coupling matrix Y„, is complex. This is required in order 
for the model to be compatible with the leptogenesis scenario of generation of the baryon 
asymmetry. In this case R can be represented as R = e* A 0, where A and O are respectively 
real antisymmetric and real orthogonal matrices. The matrix A can be parametrized by 3 
real parameters, a, b and c. We have considered the case of quasi-degenerate spectrum of light 
neutrinos, 7711,2,3 = rn u , vr? v » Am\ , Am^ , where Am\ and Anig are the neutrino mass- 
squared differences which drive the atmospheric and solar neutrino oscillations. Assuming 
that the heavy right-handed neutrinos are also quasi-degenerate in mass, Mi 2 3 — Mr, and 
that the soft SUSY breaking slepton mass terms are flavour-universal at the GUT scale, 
we have derived approximate expressions for \x — » e + 7, r — > /i + 7 and r — > e + 7 decay 
rates. Apart from the standard SUSY soft-breaking parameters (mo, do, tan/3, m^j), the 
decay rates depend on m u , Mr, on the mixing angles 812 and 823 which control the solar 
and atmospheric neutrino oscillations, on the Majorana CP-violating phases in the PMNS 
mixing matrix U and on the parameters a, b and c. We have found that for complex R, 
the branching ratios of the indicated LFV decays are considerably larger than when R is 
taken to be real: for a ~ b ~ c ~ 10 _1 , for instance, BR(fi — > e + 7) and BR(r — > e + 7) 
are enhanced approximately by a factor of 10 5 to 10 s with respect to the case of real R, 
while BR(t — > fi + 7) is enhanced by approximately four orders of magnitude. We used 
the model of leptogenesis with light quasi-degenerate neutrinos, in which the heavy RH 
Majorana neutrinos are assumed to be produced non-thermally in the inflaton decay, to get 
constraints on a, b and c. The baryon asymmetry is proportional to the product abc of the 
three parameters associated with the complexity of R. For values of the RH neutrino mass 
Mr characteristic for the leptogenesis model, the observed asymmetry can be reproduced 
for |ofec| ~ 1CT 5 . If BR(fi — > e + 7) and BR{r — > e + 7) are evaluated for values of a, b and c 
compatible with the leptogenesis constraint, the enhancement we found is approximately by 
a factor of 10 3 and 10 6 for values of sin 2 #13 = and 0.04. The corresponding enhancement 
of BR(t — > /j, + 7) is approximately by two orders of magnitude. 

Besides a, b and c, and the Majorana CP-violating phases in the PMNS matrix U, 
BR(n — » e + 7) depends also on SUSY soft breaking parameters (m , a , tan/3, myz) and 
the RH neutrino mass Mr. Given the existing experimental bound on BR(fi — > e + 7), 
our results can be used, in particular, to further constrain the space of the supersymmetric 
parameters in the case of quasi-degenerate neutrino mass spectrum. This requires a more 
detailed numerical analysis which is beyond the scope of the present work. 

If neutrinos will be proven experimentally to have a quasi-degenerate mass spectrum 
and the neutrino masses are generated via the see-saw mechanism within a SUSY GUT 
theory, the process \x — > e + 7 should be observable in the planned experiments of the next 
generation provided the supersymmetric particles have masses in the range of several hundred 
GeV. Additional constraints from data on the r — > fi + 7 and r — * e + 7 decays and from 
leptogenesis can be used to determine the matrix A and the RH neutrino mass Mr. With 
that information, the neutrino Yukawa couplings Y u could be almost fully reconstructed. 
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